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Abstract. We establish the Liouville integrability of the differential equation S{t) = 
[N, (t)], recently considered by Bloch and Iserles. Here, N is a real, fixed, skew- 

symmetric matrix and S is real symmetric. The equation is realized as a Hamiltonian 
vector field on a coadjoint orbit of a loop group, and sufficiently many commuting 
integrals are presented, together with a solution formula for their related flows in 
terms of a Riemann-Hilbert factorization problem. We also answer a question raised 
by Bloch and Iserles, by realizing the same system on a coadjoint orbit of a finite 
dimensional Lie group. 
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1. Introduction. 

In a recent paper [BI], Bloch and Iserles considered the differential equation 

S{t) = [N,S\t)], {*) 

where both S and N are real n x n matrices, S is symmetric and A is a constant 
skew-symmetric matrix. The similarity with the celebrated Toda differential equa¬ 
tion ([F], [Mo]) is hard to miss, and one might expect that both equations have 
some common properties. Indeed, from the equivalent Lax pair formulation 

S=[NS + SN, S] 

Bloch and Iserles deduced that S evolves by orthogonal conjugation: this implies 
global existence of solutions and spectrum invariance of S{t). Also, following Man- 
akov’s approach to the rigid body equation [Ma], they showed that (*) is equivalent 
to the following Lax equation with spectral parameter 

{S + zN)- = [NS + SN + zN^, S + zN], 
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from which a large collection of conserved quantities is clear, namely, the nontriv¬ 
ial coefficients in z in the expansions of tr{S -t- zN)^. Even more, they showed 
that equation (*) is a Lie-Poisson system on an appropriate phase space and, cor¬ 
roborated by numerical experiments, conjectured the complete integrability of the 
system. 

In this paper, we shall formulate equation (*) as a Hamiltonian system on a 
coadjoint orbit of a loop group, in the spirit of [DLT2]. The situation is rather 
familiar; most of the difficulty lies in guessing the loop group and a convenient 
representation of the dual of its Lie algebra. In this case, the group is LG?, the set of 
smooth loops —> GL{N, C) satisfying the reality condition g[z) = g{z) which 

admit an analytic continuation to the exterior of the unit circle with g{oo) = /, the 
identity matrix, and which are kept fixed by the involution Tj{g){z) = {g{—z)'^)~^. 
The nondegenerate pairing on arbitrary loops 

(X, Y) = <[ tr{X{z)Y{z)) 

J\z\=i 27rz ^ 

identifies the Lie algebra dual (Lg'I)* with the set of power series 5o -|- iViz , 
where the coefficients of the powers z^ ,k = 0,1,... , are alternatively symmetric 
and skew-symmetric real matrices. 

The proof of Liouville integrability of the Bloch-Iserles equation (*) takes two 
familiar forms. In one of them, embed into the larger loop group LG^ of all 
smooth invertible loops and there search for a second bracket structure, induced by 
an i?-matrix, from which the existence of commuting integrals follows. Otherwise, 
consider the Lie group LG^ x LG^, corresponding to the Lie algebra anti-direct 
sum Lq'^ © LqZ. which as a manifold can be identified with the subset of loops in 
LG^ which admit a Birkhoff factorization g = g+gZ^ with trivial diagonal part 
([PS] is a good reference for the few facts about Birkhoff factorization which will 
be used in this paper). The Lie bracket of Lq"^ 0 LqZ. is of course identical to the 
bracket defined by the ii-matrix above. Both points of view will be described in an 
intertwined fashion in Sections 2 and 3. Integrability is then proved in Section 4 by 
showing that the integrals suggested by Bloch and Iserles are indeed commuting, 
and generically span half of the dimension of the coadjoint orbit of interest. A 
formula for the commuting flows in terms of a Birkhoff factorization is the content 
of Section 5. 

The loop group scenario embeds equation (*) for nxn matrices into a coadjoint 
orbit O of generic dimension 2[n^/4] (here [x] is the greatest integer less than or 
equal to x). In Section 6, we consider a question raised by Bloch and Iserles in their 
original paper; to find a finite dimensional matrix group such that their equation 
is realized as a Hamiltonian system on the dual of its Lie algebra. In [BI], the 
authors presented an interesting algorithm to study this problem, adapted from a 
constructive proof of Ado’s theorem, which states that any finite dimensional Lie 
algebra admits a matricial faithful representation. From the loop group setup, we 
obtain an explicit realization; there is an n^-dimensional nilpotent matrix group Gf, 
faithfully realized on 3n x 3n matrices, giving rise to a coadjoint orbit isomorphic 
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to O. The prescription yielding the finite dimensional realization is simple and 
rather general. Since the orbit O is described in terms of (bi-)truncated Laurent 
expansions of its elements, one should expect that the full group is not needed to 
cover the orbit. Indeed, as we shall see, only the first three terms of a loop in LG^ 
are relevant to the computation of the coadjoint action on the orbit of interest. 
This suggests taking a quotient Gj = LG^/LG^^ by the normal subgroup LG^^ 
of loops of the form l+0{z~^): it turns out that Gf admits a simple Heisenberg-like 
representation. 

In the final section, we present an example of an infinite dimensional orbit ob¬ 
tained by the same coadjoint action. Here, the Hamiltonian which gives rise to 
equation (*) in the small orbit now induces a system of partial differential equa¬ 
tions containing some simple integral terms. 

Arieh Iserles informed us that a forthcoming paper, together with A. Bloch, T. 
Ratiu and J. Marsden, will be dedicated to a different proof of the integrability of 
the Bloch-Iserles equation, using very different techniques. 

The second author acknowledges support from CNPq and FAPERJ, Brazil. 


2. The Lie algebraic setup and the Poisson structure. 

Let G be GL{N,W) with Lie algebra 0 , and let LG be the group of loops g : 

^ GL{N,C) satisfying the reality condition g{z) = g{z). We consider the 
involutive automorphism S ; LG — > LG, given by Ti{g){z) = {g{—zY')~^ and its 
fixed point set (the stable locus) 

LG^ = {g(^ LG \ g{z){g{-z)f = /}. ( 2 . 1 ) 

Clearly, LG^ is a Lie subgroup of LG. We shall denote by (resp. LG^) the 

Lie subgroup of LG^ consisting of loops which extend analytically to the interior 
(resp. to the exterior) of the unit circle, with the additional requirement that loops 
in take the value I, the identity matrix, at oo. 

Let Lq be the Lie algebra of LG, consisting of loops X{z) = Yl°^oo^jZ^ with 
coefficients Xj G g. Then the Lie algebra Lq’^ of LG^ is the stable locus of the 
linearization a : Lq —> Lg of S, given by the formula 

a{X){z) = -{X{-z)f = ^ 9{X^){-zy (2.2) 

j 

where 9 is the Cartan involution ^ i—> Therefore, if 6 (resp. p) denote the -|-1 

(resp. —1) eigenspace of 9, consisting of skew-symmetric (resp. symmetric) matrices, 
then explicitly, 

Lq^ = { a G Lg I X 2 j G 6 , X 2 jj^i G p for all j }. (2-3) 

From the definition of (resp. LG^), it is clear that its Lie algebra Lq'^ 

(resp. Lq'^) consists of elements of the form XjZ^ (resp. Hence 

the elements in LqZ. are equal to 0 at oo and we have the splitting 


Lg"^ = LqI © LqZ. 


(2.4) 
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with associated projection maps n_|_ and n_. 

We now introduce 

G = {g & LG^ I g = g^g_^, where g- G LG^, g'_(_ G LG^ }. (2-5) 

Then an easy argument shows that for g & G, the factors g± are unique (more 
details will be provided in Section 5). Moreover, from the Birkhoff factorization 
theorem [PS], G is a dense open subset of LG^ in the natural topology. Following 
the procedure in [DLT1],[DLT2], we can endow G with a Lie group structure by 
defining the multiplication 

g*h = g+hgZ^. (2.6) 

Clearly, the map T : G —> LG^ x LG^ given by ^{g) = ^(g+gZ^) {g+,g-) 

is a Lie group isomorphism, when the image is equipped with the product group 
structure. Consequently, the pull-back of the standard Lie bracket on Lg'^ 0 LgZ 
(Lie algebra direct sum) under ^ = TeT : X i—> (Lt+X, — n_X) yields the Lie 
bracket on g = Lie{G): 

[x,y]0 = [n+x,n+y]-[n_x,n_y] (2.7) 

for all X, y G 0 . Similarly, one may use V' to obtain a formula for the adjoint action 
of G on 0 , 

Ad^{g)X = g.{I[.X)gZ^ + g+{Ii+X)g^\ (2.8) 

Remark 2.1 By standard r-matrix theory [STS], 

i? = n+ - n_ (2.9) 

is a solution of the modified Yang-Baxter equation, i.e., 

[RX,RY]- R{[RX,Y] + [X,RY]) = -[X,y] ( 2 . 10 ) 

for all X, y G L 0 '^. By an easy computation, we can show that the i?-bracket on 
Lg^ given by 

[X,Y\r = ]^{[RX,Y\ + [X,RY\) (2.11) 

coincides with [•, -jg. Thus the vector space Lg^ equipped with the Lie bracket [•, ■]r 
is identical to the Lie algebra 0 . 

To identify the duals of the Lie algebras, we introduce the following nondegen¬ 
erate invariant pairing on Lg\ 

(X,y) = / tr{X{z)Y{z))^. 

j\z\=i 27 ^^ 

= E Hx,y.,_,). 

j 


( 2 . 12 ) 
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As the reader will see, this choice of pairing is critical for what we have in mind. 
Using this pairing, we identify the algebraic dual of with the stable 

locus of the Lie algebra anti-isomorphism a* = —a: 

= { A G Lg I X 2 j+i G i, X 2 j G p for all j }. (2-13) 

Notice that the alternation of symmetric and skew-symmetric coefficients in the 
Laurent expansion still holds, now with a parity opposite to the one found in 
elements in the Lie algebra. By an easy computation making use of (2.8) and 
(2.12), we find 

Ad*Q{g)A = U_{g-^Ag+) +U+{gZ^Ag_) (2.14) 

where we have used the fact that 11^ = n_ and HI = n_|_. 

Remark 2.2 In the general context of an involutive automorphism of a finite 
dimensional Lie algebra, the author in [R] also considered the stable locus of the 
corresponding involution on the loop algebra. However, a different choice of pairing 
was used in [R], thus leading to a different identification of the dual. 

Finally, the Lie-Poisson structure for smooth functions on g* is given by the 
usual formula; 

{F^,F2}{X) = {X,[dF,{X),dF2{X)]~^). (2.15) 

3. A Hamiltonian for the Bloch-Iserles equation. 

The following result is standard. 

Proposition 3.1. For 1 < k < n, i € Z, define 

then the H^t’s Poisson commute with respect to {•,•}. The Hamiltonian equation 
of motion generated by Hki is given by 

X(z) = [u+{{Xizfz-^^+^^),X{z)] . (3.2) 

Proof. Since Ad*{X)(z) = g{z)~^X{z)g{z), the Hamiltonians are invariant under 
the coadjoint action of LG^, i.e. 

HkfiAdliX)) = Hu{X), g€LG^. 

By classical r-matrix theory, we then conclude that they Poisson commute. The 
equation of motion for is a straightforward computation. □ 

Our next task is to search for interesting finite dimensional coadjoint orbits in 
g*. For any m,n G Z_|_, define 

n 

= {Xe Lg"* I X{z) = }■ (3-3) 

j = — m 

Clearly, we can extend this definition to the case where m = 0 in which case we set 

0n =0(O,n)- 
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Proposition 3.2. The sets 0 *^ are invariant under Ad*~{g) for any g & G. 
Proof. Take X G 0 *^ For m = 0, the analyticity property of g± implies that 

Il_{g^^Xg^) = 0 , 

n — 1 

U^{gZ^Xg_)=Xr,z^ + Y,Bjz^ 

j=o 

for some matrices Bj, and we are done (we have used g-{oo) = / in deriving the 
second line above). For m G Z_|_, note that 

Ad*~{g)x = n_(5;'(n_x)5+) + n+(5z'(n+x)5-). 

Therefore, the same calculation as before applied to n_|_X shows that 

n — 1 

n+( 5 zi(n+x)< 7 _) = Xr,z^ + B,z^ 

i=o 

for some matrices Bj. On the other hand, it is easy to check that 

-1 

U_{gf\U_X)g+)=g+{0fX_mg+{0)+ Y. 

j = — m+l 

for some matrices Cj. Putting the above two expressions together, the assertion 
follows. □ 

As a consequence of this result, we obtain finite dimensional coadjoint orbits 
through the elements in 0 * and 0 *^ and the Hamiltonian equation in (3.2) re¬ 
stricts to these orbits. In the rest of the section, we shall focus on the case 0 ^, in 
which the Bloch-Iserles equation lies. 

Proposition 3.3. (a) Consider the loop Sq -|- zNq G 0 ?. The Ad*~-orbit through 

(jr 

So -|- zNq is given by the affine linear space 

Gso+zNq = {^d~,{g){So + zNq) \ g & G} = {(S'o -I- [Nq, P]) -|- zNq | P G p }. 

In particular, if Nq has simple spectrum, the tangent space of Os^+zNo is naturally 
identified with the vector space of real, symmetric matrices which are orthogonal to 
the matrix polynomials p{Nq). 

(b) The Hamiltonian equation of motion 


(S + zN)- = [U+{{S + zNf/z)),S + zN' 


(3.5) 
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generated by H 20 on the Poisson submanifold 0 ^ is equivalent to the Bloeh-Iserles 
equation 


iV = 0, 

S= [NS + SN,S]. 


(3.6) 


Proof. We have 


g-{z) = I + g'_{oo)z ^ + 0{z 


where g'_{oo) = [z = 00 ) and a simple computation obtains 
dSo+zNo = {So + [-^ 0 ) 5 -( 00 )] + zNq}. 


Since gZ^{z) = {g-{—z))'^, the matrix g'_{oo) must be symmetric and, conversely, 
any symmetric matrix P G p is g'_{oo) for an appropriate loop g-. Now, endow p 
with the usual inner product < P,Q >= tr PQ. Then the linear map Bnq ^ P ^ P 
taking P to [P, A^o] is skew symmetric and its range is orthogonal to its kernel, given 
by the set of symmetric matrices commuting with Nq. If Nq has simple spectrum, 
the description of the tangent space to the orbit then follows. The equation of 
motion for H 20 is easy to compute, once one observes that only S may vary in an 
orbit in 0 J. □ 

Remark 3.4 (a) From the proof of Proposition 3.2, it should be clear that for 
X G 0 *, the equations in (3.2) can be regarded as Hamiltonian systems on the 
coadjoint orbits of the Lie group LG^, when the dual (Lgf)* of its Lie algebra is 
identified with the collection of power series So + Niz , where the coefficients 
of the powers fc = 0, 1 ,... are alternatively symmetric and skew-symmetric real 
matrices. The advantage of using the slightly more general formulation here is that 
it puts the groups LG^ and LG^ on equal footing. Otherwise, the relevance of 
LG^ in the solution formula in Proposition 5.1 may look somewhat mysterious. 

(b) The Bloch-Iserles equation can be regarded as an isospectral deformation of 
a general n x n matrix keeping the skew-symmetric part fixed. Indeed, if we set 
2 ; = 1 in (S' -|- zN)' = [ A^S" -|- SN + zN'^, S + zN ] and let M = S + N, then a 
straightforward calculation shows that 


M = i + MM^) + 


Conversely, if M satisfies the above equation, then 


M + MM'^) + , M'^ ] 


and by a direct computation, we find (M — M^)' = 0. Now, the Bloch-Iserles 
equation preserves the spectral curve det{S + zN — w) = 0. (See Remark 4.3 below.) 
On the other hand, the Toda flow in [DLTl] also has a spectral curve, given by 
det{M-\-h{M'^ — M) — w) = 0. If we write M = S+N, then by changing the variable 
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h above in the obvious way, we can make the Toda curve coincide with det(5 + 
zN — rc) = 0. Percy Deift asked us if (3.7) might arise from a combination of flows 
associated to the Toda curve. To answer this question, note that any flow generated 
by a combination of the coefficients from the Toda curve must be of the form 
M = []![ F(M, ]. (See [DL].) Here, (is the Lie subalgebra of g consisting of 

lower triangular matrices and Ht is the projection map to ( relative to the splitting 
0 = t 0 [, and F{M, M'^) is a polynomial in M and M'^. If we denote the right 
hand side of (3.7) by X{M), then it is obvious that X{OMO'^) = 0X{M)0'^ for 
any orthogonal matrix O. However, the expression [Hi F{M, M'^), M ] clearly does 
not satisfy this invariance property. Hence the answer to the above question is in 
the negative. From a different point of view, it is also tempting to see if one can 
constrain the flows arising from the Toda curve to the submanifold Q = {M G 0 | 
M — = N } where iV Is a fixed skew-symmetric matrix. Unfortunately, the 

submanifold Q Is not a cosymplectlc submanifold in the sense of Weinstein [W]. 
(The cosymplectic submanifolds are the generalization in a Poisson context where 
we can carry out Dirac’s idea of constraining a Hamiltonian system.) Thus this 
idea also fails. 


4. Liouville integrability. 

Let symij{A, B), the ij-symmetrizer of matrices A and B, denote the sum of all 
monomials In A and B of degree i + j consisting of i A's and j B's. 

Lemma 4.1. (a) [ symij+i{A, B) , H ] 0 [ symi+ij{A, B) , B ]=0. 

(b) Let A = S be symmetric and B = N skew-symmetric. Then symij{S, N) is 
symmetric (resp. skew-symmetric) if j is even (resp. odd). 

(c) If A and B are nx n matrices, than each symmetrizer symn-e,e{A,B),i = 
0 ,... ,n, is a linear combination of symmetrizers symr,s{A,B), of smaller degree 
r -\- s <n. 

(d) For generic choices of A and B (i.e., for an open, dense set of pairs {A,B)), 
symk-e,i{A, B), k, £ = 0,... , n — 1, are linearly independent. The same is true for 
generic choices of A = S symmetric and B = N skew-symmetric. 

Proof. To simplify notation, we omit the obvious matrix dependence of symij. 
Statement (a) follows from 

symi+ij+i = A symij+i 0 B symi+ij = symij+i A 0 symi+ij B. 

The first equality states that any monomial in symi^ij^i either starts with A and 
is completed with i A's and j 0 1 B's or starts with B and is completed with 
i 0 1 A's and j B's. The second equality is obtained by taking into account the last 
matrix of each monomial. The proof of (b) is obvious. To see (c), notice that the 
claim for f = 0 (resp. i = n) is just the Cayley-Hamilton theorem for B (resp. A). 
More generally, write the Cayley-Hamilton theorem for the matrix M = A + zB, 

0 CJ 2 M "-2 + , , , + CJnM°, 
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where aj = o'j{z), a standard symmetric function of the eigenvalues of M, is a 
polynomial of degree j in z. Now collect terms in z: on the left hand side, the 
coefficient of is of the form symn-e,e] on the right hand side, it is a linear 
combination of symmetrizers of degree smaller than n. To prove (d), let be the 
matrix whose only nonzero entries, equal to one, lie along the subdiagonal of entries 
with indices (r, s) satisfying r — s = 1 and let be a diagonal matrix with entries 
along the diagonal in geometric progression c^,c^,... ,c^. Clearly, symk-e,e only 
has nonzero entries along the subdiagonal r — s = i, forming a geometric progression 
of order c^~^. Thus, to check if a linear combination of symmetrizers is independent, 
it suffices to consider the independence of subsets of symmetrizers consisting of a 
fixed index i (keep in mind that, by hypothesis, i < n), for which k — i runs from 
0 to n — £ — 1. We thus have to prove the linear independence oi n — £ vectors 
in with coordinates in geometric progressions of ratio c^,... this 

is clearly true for a generic choice of c (we thank Nicolau Saldanha for suggesting 
matrices and -B*). From analyticity, independence holds for generic choices of A 
and B, proving the first part of (d). We now show that generic independence still 
holds for symmetric (resp. skew-symmetric) choices of A (resp. B). By continuity, 
independence still holds for matrices sufficiently close to A^ and B*. In particular, 
this is still true if we keep 5* = A^ and change B* by adding a small negative 
number on its subdiagonal of entries (r, s) for which r — s = —1, giving rise to a 
matrix B^. Now, consider a diagonal matrix D so that N^, = DB^,D~^ is skew- 
symmetric. The pair = DA^D~^ and W also has independent symmetrizers of 
degree less than n and genericity for pairs {S, N) follows again by analiticity. □ 

As usual, let S and N be symmetric and skew-symmetric nxn matrices. The set 
of symmetrizers symk-e/{S, N) having degree k < n which are symmetric matrices 
(i.e., those for which j is even, by the previous lemma) splits in two sets. The first 
will give rise to the family of commuting integrals the second relates to the 
Casimirs of generic coadjoint orbits within g^: 

I = {symk-i/, 1 < k < n — 1, 0 < £ < n — 2, £ even}, 

<Sn = {symu, 0 < £ < n — 1, £ even}. 

From Lemma 4.1, for a generic choice of S and N, the map B^v is injective 
when restricted to the vector space spanned by X: in particular, the matrices 
[ symk-i/{S, N) , N ], for symk-i^i G X, are linearly independent. 

Theorem 4.2. (a) Let Nq have simple spectrum and consider an element Sq + 
zNq G g}. Then Osq+zNo is a 2[n^/4] dimensional affine linear space given by 

{{S + zN) esl \ N = No, tr{SN^) = tr{SoN^),£ even ,0 <£ <n-l} . (4.2) 


(b) The equation of motion of the Hamiltonian is given by 

S = [ symk-{e+i),e+i{S,N) , S ] = -[ symk-e,e{^y^) ^ ^ ^ = 0. 
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The Hamiltonians in the set 

{Hkt, I < k < n — 1, 0 < £ < n — 2, I even} 

are generically independent on Osq+zNo provide [n^/4] commuting integrals. 

Proof. Statement (a) follows from the characterization of the tangent space to the 
orbit given in Proposition 3.3, combined with the independence of the elements 
of C. The first formula for the equation of motion for follows by expanding 
in the equation of motion as given in Proposition 3.1. The second formula 
is a consequence of Lemma 4.1(a). Finally, by considering the second description, 
generic independence of the vector fields is a consequence of the remarks just above 
the statement of the theorem. □ 

Remark 4.3 Equivalently, we can consider the conserved quantities given by the 
coefficients of the characteristic polynomial p{z,w) = det{S + zN — wl). Notice 
that {S + zN)'^ = S — zN, so p{z, w) is an even function of z : 

N [§] 

p{z, w) = det{S + zN — wl) = EE Irk{S,N)z^'^W^-P 

r=0 fc=0 


The fact that the functions Irk Poisson commute also follows by coadjoint invari¬ 
ance. Generic independence follows from the generic independence of the Hki.. Note 
that the functions l 2 k,k, k = 1,... , [n/2] clearly depends only on N and therefore 
are trivial integrals. On the other hand, it is easy to see that {L 2 fc-i-i,fc}o<fc<[ 2 ±i]_i 

is equivalent to the coadjoint orbit invariants {tr (Hence the 
nontrivial integrals are given by {/r-fc}o<fc<[r-/ 2 ]-i,i<r-<n and this also gives a total 
of 



conserved quantities, as required. The fact that the spectral curve p{z, w) = 0 
is preserved by the Bloch-Iserles equation means that the corresponding flow is 
linearized on the corresponding Jacobian variety and that the solution can be ex¬ 
plicitly written down in terms of Riemann theta functions. We shall leave the 
details to the interested reader. (See, however, Proposition 5.2 and Remark 5.3 in 
this connection.) 


5. Solution by factorization. 

The flows associated to the integrals Hke may be described by an explicit formula. 
Let and B_ be respectively, the sets of complex numbers z together with 

z = oo for which |z| < 1, |z| = 1 and |z| > 1. Recall that ([PS]) a loop 7 G LG 
admits a Birkhoff factorization 7 (z) = 7 +(z) d{z) jZ^(z). Here, 7 + and 7 _ are 
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restrictions to of analytic functions extending to the boundary of D_|- and D_ 
taking their values on GL(N,C), the matrix d{z) is diagonal with diagonal entries 
of the form z°‘\ai G Z, and 7 _(oo) = I, the identity matrix. The only ingredient 
which is not automatic in the derivation of the formula for the flows is the proof of 
the triviality of the diagonal factor in the factorization of the loops associated to 
the initial conditions. 

As stated, the Birkhoff factorization of an invertible loop is not necessarily 
unique. What is true, and follows from an argument similar to the given in the 
proof of the proposition below, is that different factorizations have the same diago¬ 
nal factor, up to permutation of its diagonal entries. The loops we have to consider, 
however, admit a special symmetry. We follow in spirit the arguments in [GK]. 

Proposition 5.1. Let ^ be a loop satisfying the symmetry 7 ( 2 ) 7 '^ (—z) = I. Then 
its Birkhoff factorization has trivial diagonal factor, d{z) = I, and is unique. More¬ 
over, the nontrivial factors also satisfy the symmetry, 7 ±( 2 ;) 7 ^(—z) = I. 

Proof. Suppose 7 ( 2 ;) = j.^-{z)d{z)jZ^ {z), with the notation above; we have to show 
that all the integers ai are equal to 0. For a continuous function / : —> C*, 

define its winding number w{f) to be the (signed) number of turns of its image 
around the origin. The winding number is invariant under continuous deformations 
ft through functions which avoid the origin and is additive with respect to products, 
^ifg) = ^{f)Tw{g). Since 7 _|_ (resp. 7 _) extends to B_|_ (resp. ID)_), taking values 
at invertible matrices, one may deform the loop 7 to a constant loop through 
loops taking values on invertible matrices, so that rc(det 7 _(_) = rc(det 7 _) = 0. By 
additivity, rc(det 7 ) = t(;(det d) = On the other hand, since j(z)j'^(—z) = I 

and t(;(det 7 ( 2 :)) = t(;(det 7 '^(— 2 ;)), we must have t(;(det 7 ) = 0. Thus, = 0. 

Also from 'y(z)'y'^ (—z) = /, we have 

(i( 2 ;) 7 l^( 2 ;) 7 Z'^(- 2 ;)(i(- 2 ;) = 'yf\z)jf'^(-z). 

Equating diagonal entries of both sides, we obtain, for i = 1,... , n, 

{zZ\zhZ^{-z)h = 

which we denote, with the obvious attributions, by f-{z) = z~‘^°''g-^{z). Clearly, 
the hypothesis 7 _(oo) = I implies /_(oo) = 1. Suppose < 0 for some i. Then 
the function h{z), which agrees with /_ on D_ and with z~‘^°‘^ f.^-{z) on B_|_ is an 
entire, bounded function — a constant — satisfying /i(0) = 0 and h{oo) = 1. The 
upshot is that at > 0, and since X^^ at = 0, we must have at = 0 for all i. A 
similar argument involving Liouville’s theorem obtains uniqueness of factorization. 
Symmetry of 7 ± in turn follows from unique factorization applied to the equation 
Z{z) = Z+{zhL^{z) = 7i-z)~'^ = 7 +(- 2 ;)“'^ 7 _(- 2 :)'^. □ 

As we noted earlier, the vector space Lq'^ equipped with the Lie bracket in 
(2.11) coincides with g. Since the modified Yang-Baxter equation (eqn. (2.10)) is a 
factorization condition, our next result is standard ([STS]). We sketch the argument 
in order to check that the symmetries required in the previous proposition to obtain 
the Birkhoff factorization with trivial diagonal factor indeed hold. 
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Proposition 5.2. Let fke{x,z) = x^z for i even. Then the solution of the 

equation of motion 

X{z)= = [Ii+fM{X{z),z),X{z)], X(0,z) = Xo(z), 

generated by the Hamiltonian is given by 

X{t,z) = gZ\t,z)Xo{z)g_{t,z) = gf\t, z)Xo{z)g+{t, z), 
where gj^ and g- are obtained from the Birkhoff faetorization 
exp{-tfki{Xo{z),z)) = g+{t,z)gZ^{t,z). 

Proof. First, notice that, for any real t, the loop 'y(t,z) = exp{—tfkt{XQ{z),z)) 
satisfies the hypothesis of the previous proposition. Indeed, it suffices to check that, 
for even i, d{z) = —fke{Xo{z), z) satisfies the linearization 6{z)-\-6'^{—z) = 0, which 
is obvious. Thus g- and gj^ are uniquely determined from exp{—tfke{XQ{z),z)) = 
g.^.{t, z)gZ^{t, z). Differentiating the above expression, we obtain 

exp{-tfki{Xo{z),z)) fk£{Xo{z),z) 

= 9 +{t,z)gZ^{t,z) - g+{t, z)gZ^ {t, z)g-{t, z)gZ^ {t, z) 

so that 

gf^{t,z)g+{t,z) = -U^{gf^{t,z)fke{Xo{z),z)g+{t,z)) = -U^fkeiX{t,z),z) 

for X{t,z) = gf^ {t, z)XQ{z)g.^-{t, z) = gZ^ {t, z)XQ{z)g-{t, z). Taking the derivative 
of the first expression for X{t,z), we obtain 

X{t,z) = [-gf\t,z)g+{t,z),X{t,z)] = [U+fkeiX{t,z),z)),X{t,z)] . 

□ 

Remark 5.3 The Birkhoff factorization for g±{t,z) above can be solved explicitly 
in terms of Riemann theta functions. See, for example, [RS] and [DL] for details. 
Remark 5.4 In the nonperiodic Toda hierarchy, there are vector fields whose 
solution for a given initial condition Sq are especially simple at integer times. One 
example is associated to the QR factorization of exp{tf{So)) for f{x) = ln(x). 
This evolution is related to the QR iteration, a starting point of many algorithms 
to compute eigenvalues ([S],[DNT]). In a similar fashion, there are evolutions related 
to some of the Hamiltonians which are algebraically solvable at integer times, 
since the Birkhoff factorization may be performed explicitly on loops of matrices 
with polynomial entries. 
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6. A finite dimensional group for the Bloch-Iserles equation. 

We now indicate how to obtain a finite dimensional group which induces a coad- 
joint orbit diffeomorphic to Osq+zNo^ on which the Bloch-Iserles equation arises as 
a Hamiltonian system. This construction answers a question addressed in [BI], 

Represent a loop g{z) = {1+^^ + ^^ + - • •) ^ LG^ (set 5^0 = I) as the bi-infinite 
(convolution) matrix Mg which, in Fourier variables, corresponds to multiplying a 
(matrix) function defined on the circle by g{z). The matrix Mg is obtained by 
prescribing the same basis {... , z~‘^, z~^, z^, z, z^,... } on domain and range. It 
is block upper triangular, with n x n blocks indexed by {Mg)ij, equal to g-k for 
i — j = —k. Upper unipotent convolution matrices form a group Ai, and matrices 
associated to loops admitting the symmetry g(z)g'^ (—z) = I form a subgroup 
A4^. Consider the normal subgroup Alf oi on which the diagonals associated 
to fe = 1 and 2 have zero entries; such matrices correspond to loops of the form 
I + 0{z~^). The quotient A4^/A4f corresponds to ’forgetting’ diagonals associated 
to A: > 3 and is clearly isomorphic to the group Gf = LG'^/LG^^ defined in the 
Introduction. Elements in the quotient may be represented by matrices M°, on 
which diagonals for which k > 3 are set equal to zero. Now, set Ha denote the 
orthogonal projection on the three basis elements l,z,z‘^. A simple computation 
shows that is (group) isomorphic to the 3n x 3n matrices of the form 

naMgHg — this is the finite group which will induce the Bloch-Iserles equation, as 
we shall see. 

We now provide the concrete description of . Consider the group of 

real 3n x 3n matrices of the form 

fl S ^ + n\ 
g = g{S,N) = 0 / 5 , 

Vo 0 I J 


where the entries are n x n matrices. Here, matrices 5, T, 5, T are symmetric, 
N, M, A, M are skew symmetric. A simple computation gives g{S, N)g{T, M) = 
g{S + T,N + M + ■^^^), in agreement with the loop product 


,, 5 V + A T 

+...)(/ + - + 


+ ... = 


/ +-+ —^ 5 ^ + 


The Lie algebra g f and its dual g). consist of elements of the form 


/O S N\ / 0 0 0\ 

A(5,A)= 0 0 5 andA(5,A)= S 0 0 

\0 0 0 / \2A SO/ 


Here the nondegenerate pairing is (A, A) = trXA and the usual computations yield 
Ad*g^rp Up to trivial identifications, this is exactly 
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the formula for the coadjoint action of LG? through the loop S + zN described 
in Section 3. Clearly, the same Hamiltonians defined in Section 3 induce the same 
flows in this setup. 

7 . Another coadjoint orbit. 

From formula (2.18), it is clear that an nxn matrix Nq with a large kernel gives 
rise to especially small coadjoint orbits, or a concrete example, split an arbitrary 
symmetric matrix S into blocks, 

/a b u^\ 

5=16 c \ (6.1) 

\u V B j 

where a, h and c are real numbers, u and v are vectors of dimension n — 2 and B 

is a real, symmetric matrix of dimension n — 2. A direct computation shows that 
the flow 5 = [N, 5^] is equivalent to the equations 

a = 2{u, v), 
b = (v, v) — (u, u), 
c = —2(u, v) = —a, 

B =0. 

With the same notation as above, the equation 5 = [N, 5^] reads 
a = 2b^ + 2ba^ + 2b{{v, v) + {u, u)) + 2{u, Bv), 

b = —2a^ — 2ab^ - 2a{{u, u) + {v, v)) + {v, Bv) — {u, Bu), 

c = —a, (6.3) 

ii = (b^ + a^)v + bBu — aBv + {u, v)u + {v, v)v + B^v, 

V = —(6^ + a^)u — aBu — bBv — {u, u)u — {u, v)v — B^u. 

We select a simple special case. Suppose u and v are points in a possibly complex 
vector space V. Say H is a symmetric matrix (i.e., B^ = B) acting on V, so that 
its entries may be complex, and suppose that V splits in two subspaces 14 and 
Vo (we will denote them by even and odd vectors), which are real orthogonal to 
each other and are interchanged by B. In particular, we have {u, = 0, for 

all natural k, for u and v with the same parity. Finally, suppose that, at time 0, 
the entries a, b and c are 0, and the vectors u and v have the same parity. Then 
we must have that a, b and c are kept constant equal to 0 and u and v preserve 
their parity. Indeed, the requirement that a, 6 and c = —a are constant equal to 
zero is compatible with their evolutions, and the remaining equations, keeping the 
notation above to denote the real inner product on V, become 

ii = {u, v)u + (u, v)v + 

= (u, u)u — {u, v)v — B^u, 


V 


(6.4) 
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which respect parity. From uniqueness of solutions for a differential equation, what¬ 
ever solves this smaller system actually is the unique solution of the original system 
provided a, b and c stay put. 

We consider an infinite dimensional version of the above system. Take B = iD^, 
so that = —Dxx, and take u{x) and v{x) to be simultaneously either even or 
odd functions in the variable x. Let the inner product denote the usual real product 
in LF‘{dx). The evolution equations for the functions u{t,x) and v{t,x) read 

Ut = {U, V)U + {V, V)V - Vxx, 

Vt = {u, u)u - {u, v)v Uxx- 

In particular, this integro-differential system preserves the original parity and the 
reality of the initial conditions. 

Other partial differential equations which fit the format S{t) = [iV, may 
be obtained as follows. Set S = Mq the operator given by multiplication by the 
function q{x). Choose now N = DMa where D is the partial derivative 

with respect to x. Clearly, N is skew-symmetric and the evolution equation for 
the operators becomes qt = 2a{q^)x- These differential equations admit hierarchies 
of infinitely many conserved quantities, but some of these conserved quantities 
generate vector fields which are not differential equations: in a sense, the coadjoint 
orbit through the initial condition is too large. 
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